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The present qubit technology, in particular in Josephson qubits, allows an unprecedented control of 
discrete energy levels. This motivates a new study of the old pump-probe problem, where a discrete 
quantum system is driven by a strong drive and simultaneously probed by a weaker one. The strong 
drive is included by the Floquet method and the resulting quasienergy states are then studied with 
the probe. We study a qubit where the harmonic drive has a significant longitudinal component 
relative to the static equilibrium state of the qubit. Both analytical and numerical methods are used 
to solve the problem. We present calculations with realistic parameters and compare the results with 
recent experimental results. A short introduction to the Floquet method and the probe absorption is 
given. 
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I. INTRODUCTION 

The discrete energy levels of a quantum system can be 
mapped by studying the absorption from a weak harmonic 
perturbation. A coupling with a monochromatic drive 
changes the characteristics of the spectrum, an effect 
known as the dynamic Stark shift ^ . In atomic physics, 
this has a prototype in the form of an atom, driven with 
one laser, the pump, and probed with another of low 
intensity^ (see Fig. 1). Instead of the bare atomic levels, 
the probe induces transitions between the dressed states 
of the coalesced atom and pump. 

The standard case in atomic physics is a dipole tran- 
sition, where the strong drive is transverse to the static 
Hamiltonian. Similarly in nuclear magnetic resonance, 
the drive field is transverse to the static field. This limita- 
tion has been removed by the introduction of new systems 
where longitudinal drive can be included. Examples are 
Rydberg states of an atom\ Josephson qubits nitro- 
gen vacancy centers in diamond^", and semiconductor 
quantum-dots ' '^' This motivates a renewed study of 
the pump-probe physics. 

The purpose of this paper is to present calculations on 
pumped and probed qubits, where the pump has an essen- 
tial longitudinal component with respect to the equilib- 
rium level splitting. We start by a brief presentation of the 
Floquet formalism that is used to calculate the dynamic 
Stark effect of the pump field (Sec. II). The formaUsm is 
enjoying a revival due to the generation of novel systems 
allowing strong driving fields and long enough coherence 
times, e.g. superconductiirg qubits and circuits'"'^^^"^, 
and topological insulators"". The Floquet formalism" '"-'^' 
is a semiclassical method that can be understood as a 
limiting case of a full quantum mechanical picture when 
the number of quanta in the driving field is large" '. The 
resulting dressed states are called quasienergy states. We 
sketch the theory of weak probe absorption and disper- 
sion on the quasienergy states (Sec. III). This method has 
been used in the literature"^'"'''^' but, to our knowledge, 
has not been properly justified. There has been many 
experiments in the field of superconducting Josephson 



qubits"''''^''''"'"^ which could have been interpreted in terms 
of the probe absorption spectroscopy of the quasienergy 
levels. Nevertheless, only one measurement has invoked 
the method''. The general theory is applied to two level 
systems (Sec. IV). We compare numerical calculations 
with analytical approximations. We present calculations 
with parameter values that are relevant for recent exper- 
iments^' \ and compare the calculated spectra with the 
measured ones. 

The Floquet approach used here should be compared 
with an alternative method of solving the pump-probe 
problem. The conventional approach"' ' is to first find the 
steady state solution of the density matrix corresponding 
to the driven, but not probed, system. Then the probe 
absorption rate is obtained by calculating a correlator 
of the probe Hamiltonian at the probe frequency . The 
two methods are identical with the exception that the 
approximations concerning relaxation can be different. 
The Floquet approach has the benefit that the quasienergy 
structure gives additional insight and allows a simple 
analysis of the probe absorption by Fermi's golden rule. 

II. QUASIENERGY STATES 

We study a driven system described by the Hamiltonian 
H{t) 

H{t)^Ho+Hsit)- (1) 

The time-independent Hq represents the atomic system 
expressed in an atomic basis Ba = spanning the 

atomic Hilbert space Ha- The time-dependent term is 
T-periodic, Hg {t + t) = Hs (t) , and represents the strong 
driving of the atomic system. The effect of the strong driv- 
ing can be seen as a change in the energy level structure 
of the atom. By studying the absorption profile of the 
driven atom (Sec. HI), one observes that the locations of 
the spectral lines move as a function of the drive intensity. 
This dynamic Stark effect ' is depicted in Fig. 1. 
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FIG. 1. Illustration of the absorption spectrum of a strongly 
driven atom. The atom with an energy separation eq is driven 
with an angular frequency lj — 2tt/t. The drive changes 
the atomic energy (quasienergy Aq) leading to a shift of the 
resonance peak in the absorption spectrum (dynamic Stark 
shift). 



The time-dependent Schrodinger equation is 
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Due to the periodicity of H{t), this is now analogous 
to the one- dimensional Bloch's problem of solid state 
physics^''. Within the analogy, the solution of Eq. (2) can 
be expressed in the form 



u{t)) 



(3) 



where the state \u{t)) is r-periodic and e is called the 
quasienergy-' ' . 

In order to solve Eq. (2) by using (3), the Floquet 
method takes advantage of the periodicity of H{t) and 
\u{t)). The atomic Hilbert space Ha is expanded with 
T-periodic functions, Hr, spanned by the basis Br = 
{\n) ,n € W,: {t\n) = exp{inLdt)}, where w = 2tt/t. This 
composite Hilbert space Ha Hr is referred to as the 
Sambe space"'. This expansion allows the representation 
of the periodic quantities in terms of time-independent 
coefficients 
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Eqs. (4) and (5) can also be seen as the Fourier series 
representation oi Hfj„i{t) and \ua{t)). 

In the Sambe space, the Schrodinger equation (2) re- 
duces to 
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which can be understood as a time-independent (Floquet 
matrix) eigenvalue problem 



Hp \u) = e \u) 



(7) 



with an infinite rank. The quasienergy e can now be under- 
stood as the combined energy of the atom and the driving 
field, and it is analogous to the quasimomentum in solid 

state physics ^ Quasienergy states \u) = X^o-n ^^^"^ l"';'^) 
are obtained as the eigenvectors of the Hermitian matrix 
Hp, and thus form a complete and orthonormal basis of 
the Sambe space. Generally, the eigenvalue equation (7) 
has to be solved numerically by using an appropriate 
truncation. In some limits approximate analytic solutions 
can be found (see Sec. IVA-IVB). It is worthwhile to 
notice the beauty of the Floquet approach, it reduces 
time-dependent problems to static ones, that is, to time- 
independent eigenvalue problems. 

The quasienergies obtained from the eiqenvalue equa- 
tion (7) have a periodic structure. Corresponding to a 
quasienergy e^, there is an infinite set of solutions with 
quasienergies er.n = + nHoj, where n is an integer. The 
corresponding eigenstates | trivially obtained 

from each other as they produce the same \^{t)) in Eq. 
(3). Therefore it is sufficient to solve numerically the 
states in a single quasienergy interval of width hiu, which 
is referred to as a single Brillouin zone. The number of 
such states equals the number of basis states of the atomic 
Hamiltonian Hq. Physically the quasienergy states can 
be interpreted as the atomic states being entangled with 
the driving field containing different number of quanta. 



III. PROBE SPECTROSCOPY OF 
QUASIENERGY LEVELS 

The quasienergy spectrum can be studied in terms of 
absorption from a weak perturbation, similar to the time- 
independent quantum systems. As a demonstration of 
the power of the Floquet method, we derive the transition 
rate between two quasienergy states in a similar fashion 
to Fermi's golden rule. However, we generalize the time- 
independent results, typically derived by using harmonic 
perturbation, by allowing the probe Hamiltonian to be 
quasiperiodic 



Hv{t) = Fp{t)Fs{t) + Fp{t)Fs{t) 



t 



(8) 



(Hp)<"r' 



Here, Fp is rp-periodic and Fa, is r-periodic. When 
the periods are incommensurate, the probe Hamiltonian 
Hp(t) is not periodic, in spite of consisting of products 
of periodic operators. The quasiperiodic form (8) allows 
various realizations of the probe' We take the rp- 
periodic part of the probe (8) to have the harmonic form 
Fp{t) = Ape~"^p*, with the amplitude ^p and the angu- 
lar frequency wp = 27r/rp. More general functions can 
be decomposed into Fourier series, where each term can 
be treated independently of the others. The r-periodic 
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part of the probe becomes time-independent in the Sambe 
space, and can be represented in the quasienergy basis as 



(9) 



where the summations go over aU quasienergy states (aU 
Brillouin zones). The matrix elements 



pq 



Up{t) 



Fs{t) Ug{t)) = 



(10) 



do not depend on time and are easy to implement after the 
numerical solution of the Floquet eigenvalue problem (7). 



A. Golden rule for transitions between quasienergy 

state 

We include the probe Hamiltonian (8) into the 
Schrodinger equation (2), and look for a solution in the 
form 
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where the summation goes over all quasienergy states. 
This leads into a differential equation for the probability 
amplitude a f (t) to be in the state / 
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where 00 fi = (e/ — ei)/ti denotes the transition frequency 
between the quasienergies. 

We now assume that the probe amplitude is small. 
When the matrix elements Ffi are small in comparison 
to the characteristic energy quantum hiop, it is sufficient 
to make a perturbative expansion in Eq. (12) up to first 
order in Ap '^. The integration over the broadened final 
state leads to a finite steady state rate Vi^f for the 
transition from the initial \^>i(t)) to the final quasienergy 
state |^/(t)). We assume that the transition rate is small 
in comparison to the broadening of the final state. The 
summation over all final quasienergy states gives the 
absorptive {ujfi > 0) transition rate 
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where 7/^ = 7i + 7/ is the sum of the widths of the ini- 
tial and finals states, which both are assumed to have 
Lorentzian form. We have also summed over all initial 
quasienergy states weighted by their steady state occu- 
pation probabilities pj. We have neglected all but the 
resonant term in Eq. (12). This is analogous to the 



rotating wave approximation (RWA) where the rapidly 
oscillating terms are assumed to average out in the steady 
state dynamics. 

The result (13) can be named as Fermi's golden rule for 
transitions between quasienergy states, as it is analogous 
to the result obtained between energy eigenstates " . The 
transitions between the quasienergy states occur when the 
corresponding quasienergy difference equals the energy 
quantum of the probe: e/ — = hujp. The magnitude 
of the transition is proportional to the squared matrix 
element l^/ip. It is worthwhile to note that the transi- 
tion does not occur between the atomic states \a), but 
between the quasienergy states \u). As seen by the atomic 
system, multiple strong driving quanta can participate 
in the process since the quasienergy states can be in any 
Brillouin zone: e/o ~ ^io + {1^ ~ n)huj = hup. Yet, within 
the first order approximation and with harmonic pertur- 
bation, only one probe quantum can be exchanged in the 
transition process. 



B. Relation to the spectrum of the probe field 

The transition rate P can be expressed alternatively 
using correlation functions. By applying the completeness 
of the quasienergy states ju^) one finds that V (13) is 
equal to the absorption spectrum 



\At 
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(14) 



at the probe frequency. The same expression is also ob- 
tained by formulating the probe absorption spectroscopy 
using the input-output- formalism ' In the correspond- 
ing emission spectrum, the order of the operators in the 
correlator is interchanged. The correlator approach (14), 
usually calculated using numerical integration of the mas- 
ter equation, gives the same information about the loca- 
tions and widths of the spectrum peaks as the Floquet 
approach (13). Nevertheless, the possible resonance shifts 
(i.e.. Stark and Bloch-Siegert ' ' shifts) or the magnitudes 
of the resonances are cleanly explained in the Floquet 
method with the quasienergy structure'' giving physical 
insight on the composition of the driving field and the 
system. 

In the linear response theory, the absorption spectrum 
is given by the imaginary part of the generalized probe 
susceptibility a{LOp) — a'{LUp) + ia"(a;p)'^-, a function 
which determines the dynamics of the system under per- 
turbation. This is generally referred to as the fluctuation- 
dissipation theorem. Thus, the imaginary part a" (up) be- 
comes directly proportional to the absorption V{ujp) (13) 
of the system. Moreover, the real part a'{ujp) (disper- 
sion) , which makes the phase shift of the response, can be 
obtained analytically from a"{u!p) [in practice Eq. (13)] 
by using the Kramers-Kronig -relations. This way one can 
solve both the absorption and the dispersion by exploiting 
the quasienergy structure. 
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C. Extensions 

The golden rule (13) is a perturbative result in the 
perturbation parameter A = FfiAp/hujp. The second 
order contributions become significant when A is com- 
parable to unity. With such large transition strengths, 
the original quasienergy structure becomes altered by the 
probe field. Instead of calculating the higher order expan- 
sions in A, we propose the use of the generalized Floquet 
method ' .It allows the calculation of the quasienergies 
of a Hamiltonian having two, or more, driving fields with 
arbitrary driving amplitudes Aj and (incommensurate) 
frequencies ujj. The resulting quasienergy structure is 
'quasiperiodic', which in the case of two driving fields 
means that £r,n,m — Sr + nhui + 7nhLj2. 

The detailed analysis of the two-mode quasienergy 
structure provides a quantitative method to study, among 
others, the validity limit of the first order expansion lead- 
ing to the golden rule (13). A comparison can be made by 
studying the differences between the quasienergy levels 
calculated with and without the probe field. One can 
say that the golden rule consideration is not valid if the 
results differ outside the expected locations of the weak 
probe resonances (anti-crossings), or if these locations are 
shifted. Details of the generalized Floquet method ap- 
plied to the strongly driven and weakly probed quantum 
two-level system are given in Appendix. 
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FIG. 2. Schematics of the longitudinally driven and weakly 
probed two-level system both in the temporal space and the 
quasienergy space. The temporal space; A two-level system, 
whose energy splitting e{t) oscillates sinusoidally around the 
mean eo with the period r. The quasienergy space; Trans- 
formation to the Floquet formalism results in the temporally 
static quasienergy levels that repeat in energy with the period 
hu). The weak probe field (blue arrows) acts between the 
atomic states or between the quasienergy states. In the case 
of Up < u, the probe resonance is met when Aq = fiwp or 



■ hu-p. 



IV. APPLICATION TO A TWO-LEVEL 
SYSTEM 

We give an example on the probe spectroscopy of 
quasienergy states by studying a two-level system un- 
der a strong longitudinal drive ''''^''^'-^ (see Fig. 2). In 
similar systems^'^'^"''^ one has previously considered 
the rotating wave approximation (RWA) and the Landau- 
Zener-Stiickelberg (LZS) approach ' ', whose point of view 
is in the discretized, 'stroboscopic', evolution of the pe- 
riodically driven qubit in the temporal space. In the 
LZS-approach, the inclusion of an additional probe field 
is complicated. In contrast, we concentrate on the pos- 
sibility to directly map the quasienergies by studying 
absorption from the probe. 

We assume the Hamiltonian 

TT/ \ A~ A , Ap 

= -T^Oy. + -T^Ox + COS(wi)CT3 + — C0S(wpi)cr2, 

(15) 

where the operators ^x,y,z denote the Pauli spin matrices. 
Here the first two terms form the atomic part Hq, which 
consists of a static energy splitting Eq and a tunneling 
amplitude A. The third term is the strong drive Hs{t)- 
Together with the first term, this implies that the level 
spacing e{t) = eo+Acos{ujt) (neglecting A) oscillates with 
amplitude A and frequency uj = 2'k/t. After transforming 
to the Floquet formalism, this is reflected in the periodic- 
ity in the quasienergy, see Fig. 2. In a two-level system. 



we define the quasienergy splitting Aq — £-)_ — £_ as the en- 
ergy difference between the two quasienergy levels within 
a Brillouin zone. This, together with the /zw-periodicity, 
includes all relevant information about the energy level 
structure of the driven two-level system. The fourth term 
in (15) is the probe Hamiltonian Hp{t). It is assumed to 
act in the same direction as the drive, but with a small 
amplitude Ap and a different frequency wp. The same 
direction can be arranged, e.g., by coupling the probe to 
the system through the same channel as the strong drive. 
For simplicity, we consider here a purely rp-periodic probe 
Hamiltonian (8) by setting Fs{t) = az- Rcf. gives an 
example of the probe absorption spectroscopy in the case 
of a non-trivial quasiperiodic probe. 



A. Choice of basis 

As was discussed in Sec. II, the infinite (in rank) Floquet 
Hamiltonian has to be truncated before its eigenproblem 
can be solved. The accuracy of the truncation is depen- 
dent on the choice of the atomic basis Bp^. In the case of 
a strongly driven two-level system, there are two natural 
choices for the basis, the adiabatic and the diabatic bases 
(see also Ref. 40). Here, the eigenbasis of CTz in (15) is 
called the diabatic basis. It holds the implicit assumption 
that the tunneling amplitude A is a small perturbation, 
A/fiw ^ 1. Another choice for the basis is the eigenstates 
of the static Hamiltonian /fp- This is referred to as the 
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adiabatic basis, which works the best when the tunnehng 
amphtude A is not just a smah perturbation, but of the 
same order as huj and Eq. 

In the presence of substantial driving, one way to 
decide the basis preferable for the calculations is to 
study the LZS-dynamics'"''*'* of the driven qubit H{t) — 
Ho + Hs{t). The probabihty of Landau-Zener (LZ) tun- 
nehng between the adiabatic eigenstates is given by 

Plz = cxp (^-2TrA'^/4:hujy/A'^ - e^j for A > Eq, other- 
wise Plz is small. If Plz is small, the adiabatic basis is 
the natural choice for the basis in quasienergy calculations. 
In the opposite case where Plz ^ 1, the diabatic basis 
states are closer to the eigenstates of the Floquet Hamil- 
tonian, and thus appropriate for quasienergy calculations. 

In the following analytic calculation of the quasienergy 
structures, we use the adiabatic basis when A < Sq and 
the diabatic basis otherwise. After solving the quasiener- 
gies, we consider the probe induced transitions between 
quasienergy states in the diabatic basis using the RWA. 
It is important to note that whereas the approximate 
results are basis dependent, all exact results (such as the 
numerical quasienergies) are not. Nevertheless, the size of 
the truncated Floquet Hamiltonian required for accurate 
results can have strong dependence on the chosen atomic 
basis. 



B. Quasienergy states 

We neglect the probe and dissipation, and consider 
only the strongly driven qubit H{t) = Hq + Hs{t). We 
do a transformation into a non-uniformly rotating frame 
with H' ^ U^HU + \h{dtU^)U, where U{t) is a unitary 
time-dependent rotation 



lJ{t) = exp 



-i^ sH-t)a.^ 



(16) 



The operation removes the strong drive in the z direction 
at the expense of generating in the x direction all harmon- 
ics nuj with relative weights A J„(^//iw)/2. According 
to Sec. II, all T-periodic entities are time-independent 
in the Sambe space and can be expressed in the matrix 
notation ^ . 

A resonance between the strong drive and the qubit 
is seen in the Sambe space as a pair of states that are 
nearly degenerate. Here we assume that the contribution 
of the non-resonant states to the resonant coupling is 
small. Thus, one can rely on the RWA and ignore all but 
the resonant states and the direct coupling between them. 
We choose one pair of the resonant states, resulting in 



H^WA 
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(17) 



There is an infinite amount of other similar pairs that 
are just copies of Eq. (17) shifted in energy due to the 
periodicity of the Floquet matrix Hp. 



The diagonalization of -ffRWA produces the quasienergy 
difference A^^A 



Hrwa 



A, RWA 

2 'I 



^^/(eo-nM^ + AVM^). 



(18) 

The diabatic-basis RWA is accurate if A/Zioj ^ 1. By 
following the generalized van Vleck perturbation the- 
ory^'''-"'^' , the RWA result (18) can be corrected with 
higher-order terms in the perturbation parameter A/huj. 
The second order correction^ affects the locations of 
the strong driving resonances: Eq = ntvxj — 5. We call it 
the 5 shift. In the diabatic basis, the explicit expression 
for 5 shift is 



oo 
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(19) 



The corrected quasienergy splitting is then 



^RWA+i5 



nfiw)2-fA2j2 ( AY 



(20) 



In the diabatic basis, the shift is the most important 
at small amplitude A and at small n. This implies that 
the i5d-shift vanishes with moderate driving amplitudes' ' , 
that is, when the diabatic basis is the most natural choice 
for the basis. 

In the adiabatic basis, one first diagonalizes Hq and 
then transforms H(t) to the non-uniformly rotating frame 
with a time-dependent transformation analogous to (16). 
The resulting Floquet matrix has exactly the same struc- 
ture as in the diabatic basis, but the diabatic diagonal 
energy Eq is replaced by hujQ = Y^e§ -t- A^ and diabatic 
coupling strength 



AJniA/huj)/2 



nhujA 
2eo 



Jn 



A eg 



(21) 



The adiabatic resonance condition can then be written as 
hujo = nhuj — (5a, where the adiabatic (5a shift is calculated 
with the formula (19), but using the adiabatic coupling 
strengths (21) and the diagonal energies fiujQ. 

In Fig. 3, we have shown the comparison of the numer- 
ically and analytically calculated quasienergy landscapes 
in the ea — A plane. The adiabatic basis (black dashed) is 
applied when A < £q and the diabatic basis (red dashed) 
otherwise. The analytic quasienergies agree well with the 
corresponding numerical ones when the effects of tunnel 
amplitude can be handled with the perturbation theory, 
cf. Fig. 3(a)-(b). But, the generalized van Vleck pertur- 
bation theory becomes insufficient"" if the fraction A/ A 
becomes large enough, and simultaneously A/hui > 1. In 
this limit, the calculation of the quasienergies is necessar- 
ily numerical. The breakdown of the analytical approach 
is demonstrated in Fig. 3(c)-(d). 
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FIG. 3. Landscapes of the quasienergy Aq in the eo — A plane with diflFerent values for the tunnel amplitude: (a) A/huj = 0.10; 
(b) A/haj = 0.37 corresponding to an experimental realization ; (c) A/hu = 0.84 corresponding to another experimental 
realization ; (d) A/tku — 1.50. The energy scale is the same in all panels. For the contour lines, we use the following color 
coding: numerical (solid gray), analytic in the diabatic basis Eq. (20) (dashed red), and analytic in adiabatic basis, that is, 
the adiabatic version of Eq. (20) (dashed black). The numerical contour lines are spaced by 0.10 ficj, except for two contours 
in the panel (b), where the comparison between the numerical and analytical contour lines is done with Aq/ficj = 0.092 and 
Aq/fioj = 0.918 since they correspond to the values of the experimentally measured weak probe resonances . 



C. Weak probe transitions 

We now discuss the probe resonance condition and 
the probe transition elements Ffi (10) in terms of the 
diabatic basis and the RWA. This kind of treatment is 
adequate for the essential physical insight. We follow 
the same procedure as in calculating quasienergies. The 
transformation U{t) (16) does not change the probe part 
of Hamiltonian (15). In the subsequent transformation to 
the Sambe space, the r-periodic part of the weak probe 
obtains the form Fs = CTz €51, where I denotes the infinite- 
dimensional identity matrix. 

As the strongly driven part of the total Hamiltonian is 
truncated into a two- level system (17), it is reasonable to 



make the same reduction for the weak probe part. The r- 
periodic part of the weak probe becomes simply Fg — cr^, 
operating between the resonant basis states of the RWA 
Hamiltonian (17). In the diagonalization of the strongly 
driven part of the Hamiltonian, the perturbation matrix 
gets a non-diagonal form 



S — A RWA 
9 



AJ„ 



'A_^ 

^ huj J 

^RWA 



(22) 



expressed directly in the basis of quasienergy states \u+) 
and with the energy splitting A^^^ (18). The lon- 
gitudinal weak probe itself would not induce transitions 
between the non-driven diabatic eigenstates, but the ro- 
tation to qubit eigenbasis (cx A) and the dressing of the 
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FIG. 4. Comparison between the numerical (solid line) and 
RWA analytical (dotted line) transition amplitudes l-F/iP (24), 
calculated with the parameters of Fig. 3(b). The transition 
amplitudes are picked by following the corresponding resonance 
conditions (23), i.e., from the parametrized line eQ{A) where 
Aq = 0.092/lw or Aq = 0.918/lai. The curves are labeled by 
the index n in Eq. (17). 



strong drive [oc Jn{A/huj)] have such a effect that probe 
transitions become possible. This is formally seen as the 
non-zero transverse cTx term in Eq. (22). 

In the general case, the resonance condition for the 
probe transition is 



hujp — kfujj. 

{k + l)huj — hojp, 



(23) 



where fc = 0, 1, 2, 



The quasienergy difference Aq 



is chosen so that kuj < uip < (fc-l-l)cL'. 

is defined as the difference 
between two consecutive quasienergy levels. We consider 
now the case k = (cf. Fig 2). Thus, the weak probe 
transitions are possible when the quasienergy difference 
Aq = Hujp or Aq = huj — fiwp, and the corresponding 
transition matrix element is non-zero. In Fig. 3, the 
resonance condition is shown as highlighted contour lines 
(dashed lines) . If the tunneling amplitude A is so large 
that the minimum quasienergy difference |AJ„ (A/hu) \ is 
larger than the probe energy hujp, there are no resonances. 
In Fig. 3(a), the tunneling amplitude A is small and the 
resonances are continuous lines in vertical direction, but 
as the value of A is increased the resonances curve and 
close, cf. Fig. 3(b)-(d). 

The matrix element (10) for the weak probe transition 
between the quasienergy states and is directly 
the non-diagonal element in (22) 



(24) 



We are interested in the transition element when 
the weak probe is (nearly) resonant, that is, Aq ^ 



hujp 



Thus, the transition amplitude depends only on the cou- 
pling strength AJ„{A/hLL)) of the two uncoupled energy 



levels in Eq. (17). The comparison between the numeri- 
cal (solid) and analytical (dotted) transition amplitudes 
is shown in Fig. 4. It is calculated by following the 
weak probe resonances (23). The agreement between the 
numerical and analytical results is good by taking into ac- 
count that the chosen parameters are close to the validity 
boundary of the RWA. 

To calculate the transition rate V (13), in addition 
to the quasienergies and the quasienergy states, one 
needs the dephasing rate = 7 and the populations 
Pi of the quasienergy states. We estimate them by fol- 
lowing Refs. 2(J and 42 that apply the Floquet-Born- 
Markov-formalism-"'^^'"'""'"^, which successfully merges 
the Floquet method and detailed coupling to the en- 
vironment. First, one constructs the master equation 
for the strongly driven qubit coupled to the environ- 
ment through the az operator, i.e., via the matrix el- 
ements Xai3n = (wq,o lo'z I| U/3,n)- The quasicuergy 
states \up^n) are employed to calculate the above matrix 
elements Xafsn- Here this is done numerically, but it can 
also be done analytically with the RWA or with the second 
order Van Vleck-correction, within their validity ranges"'". 
The environment is modeled with a continuum of har- 
monic oscillators, i.e., a thermal bath characterized with 
Nai3n = \Gai3n {coth [h(ea — ^0 + nuj)Hu! /2kT] — 1} and 
the Ohmic spectral density Gapn = G{ea — + nuj) = 
K{ea —5/3 + nuj). Finally, the coefficients in the master 
equation are averaged over the period of the strong drive, 
in order to bring them into time-independent form'^'' *^'' 
(secular approximation, moderate rotating wave approxi- 
mation). The result is analytically solvable in the steady 
state limit '", from which the dephasing rate 7 and the 
population p_ (p_|_ = 1 — p_ ) are derived 



CO 

7 = 7r ^ (2iV_+„ + G_+„)X2^„ + 4A^__„X2 

n— — 00 



p- 



(25) 
(26) 



The numerically calculated transition rates V (13) are 
shown in Fig. 5 in the Sq — A plane. 

The total line-shape (13) encodes the information on the 
quasienergy structure at the locations of the resonances 
and on the transition amplitudes in the magnitudes of 
the resonances. By comparing the line-shapes in Fig. 5 
with the quasienergy structure of Fig. 3, one observes the 
faithful mapping of the energy landscape. The maximum 
value of the transition element l-F/iP (24) depends on 
the tunneling amplitude A. If A is large enough, the 
maximum is reached when AJn{A/huj) — hujp. The 
transition element cannot obtain larger values since then 
the resonance condition is not anymore valid, see Eq. (18) 
and Fig. 5(b). With smaller tunneling amplitude A, the 
maximum of the is directly set by the maximum of 

Jn{A/tujj), see Fig. 5(a). The weak probe signal vanishes 
at the zeros of the Jn{A/fkd), which are related to the 
coherent destruction of tunneling . This is seen in Fig. 5 
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1 2 EQlhuj 3 1 2 solhx) ^ 

FIG. 5. The transition rate V (13) of the strongly driven and weakly probed qubit presented as a gray-scale plot in the eo — A 
plane. The parameters are the same as in the corresponding panels in Fig. 3. In addition to those, the parameter k describing the 
Ohmic spectral density is chosen so that "//uu = 0.016 in the absence of driving and detuning in Eq. (25) and /3 = hui/kT = 2.24. 
The discontinuities of the lines are a consequence of the roots of Jn{A/hu) related to the coherent destruction of the tunneling. 
The V versus A plot along the vertical dashed line in panel (b) is shown in Fig. 6. 



as discontinuous resonance lines, although the underlying 
quasienergy resonance conditions are continuous lines (a) 
or closed curves (b)-(d). 



D. Relation to the spectrum of the probe field 

In the case of the two-level system (15), the spectrum 
as a function of the correlator (14) takes the form 

S{UJP) = ^ {^.{t)aM) e"^^*dt . (27) 

Noteworthily, this spectrum is not the one commonly 
calculated from the transverse correlator ((T_(r)(T+(0)), 



natural to the atomic systems coupling to the environ- 
ment through the (transverse) dipole moment. In Fig. 6, 
we have compared the line-shape V calculated by using 
numerically implemented Floquet method (solid line) and 
weak probe response (circles) S{uj-p) (27), obtained by 
solving the steady state master equation. The correlation 
function approach (circles) agrees very well with the tran- 
sition rate calculated with the numerical Floquet method 
(solid), which further validates the method of the probe 
spectroscopy of quasienergies. The slight differences be- 
tween the two methods can be traced back to the different 
approximations concerning relaxation and dephasing. In 
contrast to the detailed Floquet-Born-Markov-formalism, 
the master equation of the qubit corresponding to (27) 
includes simply the standard relaxation and dephasing, 
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4 A/hcj 6 



FIG. 6. The transition rate V of the strongly driven and 
weakly probed qubit and the corresponding spectrum S(cjp), 
both scaled with their maximum value. The parameters are 
eo/hu) = 1.05, A/fiij = 0.37,and cjp/oj = 0.092. The pa- 
rameter K describing the Ohmic spectral density is chosen so 
that 7/0; = 0.016 in the absence of driving and detuning in 
Eq. (25) and P = hw/kT = 2.24. The transition rate V is 
calculated with the numerical Floquet method [solid line, ver- 
tical projection from Fig. 5(b)], which is contrasted with the 
spectrum S{u)p) (27) (circles) of the driven qubit (I/T2 = 7 
and l/Ti = 7/2) at the weak probe frequency. 




FIG. 7. The probe absorption V calculated as a function 
level spacing eo and driving amplitude A. The parameters 
are u/2-k = 7.0 GHz, A/fio; = 0.37, a;p/cj = 0.092, and 
T = 150 mK. The parameter k describing the Ohmic spectral 
density is chosen so that 7/0; = 0.045 in the absence of driving 
and detuning in Eq. (25), corresponding the experimental 
estimate for qubit dephasing. This plot should be compared 
with the experimental plot in Ref. 8. For the comparison we 
have given the axis scales also using the units of this reference. 



with rates 1/Ti and I/T2, respectively. 



E. Comparison with experiments 

We have also interpreted two recent experiments in 
terms of probe absorption of quasienergy states. The 
experiment by Wilson et al. uses a Cooper-pair box and 
the experiment by Izmalkov et al.'' uses a flux qubit, but 
both can be described by the Hamiltonian in Eq. (15). 

Figure 7 shows the calculated probe absorption cor- 
responding to the experiment of Wilson et al . The 
parameters are the same as given in Ref. except that 
we have not included the extra broadening caused by 
low- frequency fluctuations in the gate charge rig. We 
have used the same parameters also in Figs. 3(b) and 
5(b), except that the line- width 7 is almost three times 
larger than in Fig. 5(b). The resonances in Fig. 7 still 
have the same characteristic features as in Fig. 5(b), but 
they are not as clear because of the larger line width. 
Figure 7 should be compared with the experimental plot 
in Ref. which, however, has the extra broadening that 
wipes out some of the features. In the same reference, the 
experimental data was successfully compared with theory 
by using RWA, which is still sufScient at the parameter 
values of the experiment (see Fig. 3). 

Figure 8 shows the probe absorption calculated with the 
parameters corresponding to the experiment of Izmalkov 
et al''. The same parameters are also used in Figs. 3(c) 
and 5(c), except that the probe frequency wp/w = 0.005 



is much smaller than in Fig. 5(c). Now, the elliptical 
shape of the resonances is not resolved because of line 
broadening, but the discontinuities of the resonances re- 
main. The 5d shift (19), which is a signature of the RWA 
breakdown, is clearly visible as the bending of the res- 
onances as a function of driving amplitude A. The i5d 
shift is enhanced near — and at small n^-'\ The plot 
should be compared with the experimental plot in Fig. 3 
by Izmalkov et al. taking into account that it is a phase 
plot instead of an absorption plot. Both plots reveal the 
same quasienergy landscape, as the resonances are visible 
as bluish lines and the discontinuities as yellow crosses 
in the phase plot. In the same reference the experimen- 
tal results are interpreted as LZS-interferometry " which 
produces oscillations of the qubit population. 

The discussed experiments reveal information about 
quasienergy landscape, but suffer from noise that prevents 
the observation of individual contour lines. We point out 
Ref. f> as an example of an experiment where individual 
contour lines are clearly seen. Another difference in this 
experiment is that the modulation of the energy is non- 
sinusoidal, leading into genuinely quasiperiodic probe, 
in contrast to Eq. (15). The Floquet analysis at the 
parameters of this experiment was reported in conjunction 
with the measurement (see Supplementary Information 
of Ref. G). 
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FIG. 8. The probe absorption V calculated as a function 
level spacing eo and driving amplitude A. The parameters 
are uj/2-k = 4.15 GHz, A/fej = 0.84, uj-p/i^ = 0.005, and 
T = 70 mK. The parameter k describing the Ohmic spectral 
density is chosen so that 'y/uj = 0.17 in the absence of driving 
and detuning in Eq. (25), corresponding to the experimental 
estimate for qubit dephasing. This plot should be compared 
with the experimental plot in Ref. . For the comparison we 
have given the axis scales also using the units of this reference. 



V. CONCLUSIONS 

We have presented a method to map the quasienergies 
of a driven quantum system by using a weak probe. We 
made the derivation with a general form of the probe 
Hamiltonian, but applied it to simple cases in order to 
gain physical insight. Provided that the quasienergy exci- 
tation has a long enough life time, the spectroscopy en- 
ables an accurate mapping of the quasienergy structures". 
The results rely on first order perturbation expansion 
in the probe amplitude. We also suggested the general- 
ized Floquet method as a possible way to go beyond the 
perturbative-probe approximation. 

The detailed discussion about the strongly driven and 
weakly probed qubit shows that, with certain parameter 
values, analytical results may be obtained for the weak 
probe resonances and the transition amplitudes, thus re- 
sulting both the absorption and dispersion of the probe 
response, i.e., the generalized probe susceptibility. Oth- 
erwise, numerical calculations are a necessity. However, 
relying only on proper matrix truncation and inversion, 
the solutions are numerically stable and simple to find. 
We noted that the accuracy of the analytic, and to some 
extent the numerical, calculation is dependent on the 
choice of the atomic basis. Indeed, the detailed study 
of the transition from the adiabatic to the diabatic be- 
haviour would be interesting and possible by using the 
probe absorption spectroscopy of quasienergies. 

We reinterpreted two recent experiments'^' '. Although 
the estimated life-time of the quasienergy excitations in 



the referred experiments were too short to distinguish the 
quasienergy contours, we were able to point out features 
in the measured responses that stem from the underlying 
quasienergy landscape. 
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Appendix: The generalized Floquet method 

The generalization of the Floquet method was devel- 
oped in Ref. 3G. It enables the handling of bi- or poly- 
chromatic driving fields in a way similar to the monochro- 
matic case. Here, we use the two-mode Floquet method 
in the analysis of the strongly driven and weakly probed 
qubit. We assume the bichromatic Hamiltonian defined in 
Eq. (15). In the generalized Floquet picture, the solution 
of the time-dependent Schrodinger equation is given in 
the form 



-\ttjh 



u{t)) 



-ih—+Hit)]\u{t))=e\uit)) 



(A.l) 
(A.2) 



where the quasienergy state \u{t)) is also bichromatic 
and the quasienergies are quasiperiodic e — er,ni,n2 — 
Cr + hniuj -\- hn2Ujp. 

To take advantage of the periodicity, we express the 
Hamiltonian (15) and the state \u{t)) using a 'double' 
Fourier series representation: 



H{t) 



ni— — oo cr,a' 
712 — — oo 



^ E ^e'("i'^+"^'^p)*c(,"i)'("^) |a) . (A.4) 



ni— — 00 cr 
712— — 00 



Similar to the case of the single-mode Floquet method [see 
Eq. (7)], we get a time- independent eigenvalue equation 



Hf, \u) = e \u) 



(A.5) 



The Hamiltonian (15) can be expressed in terms of the 
sub matrices = i (eoo-^ -I- Ad-^;), H[='=i1 = ^ct^, and 

Hit) = H["1-HH[±i1 (e'"* e"'"*) 

_^ B[±i,o] ^gi^pt _^g-ic.pt^ ^ (A.6) 
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FIG. 9. Quasienergy in the generalized Floquet method, (a) Quasienergy landscape in the eg — A plane with the parameters 
A/hu} = 0.37, iop/uj — 0.10, and Ap/fko = 0.20; The gray line shows the weak probe resonance condition (23) deduced from the 
corresponding single-mode quasienergy landscape [Fig. 3(b)]. (b) Projection of (a) at A/fuo = 5.6, denoted with arrows and 
dashed line. The quasienergy er,n,m (black) is periodic so that er (blue) is shifted by nhoj + mhojp, where n,m — Q, ±1, ±2, . . .. 
The red dash-dotted line shows the corresponding quasienergy calculated with the single-mode Floquet-method. (c) Comparison 
of the two-mode quasienergies with the probe amplitude Ap/huj — 0.20 (solid) and 0.40 (dashed). Magnified view of the box in 
panel (b). 



The two-mode Floquet matrix Hp^ of the Hamiltonian is 
given as an infinite dimensional matrix ''' 



Hp - Ihujp Bill 
Bl-i] Hf BW 
Bl-i] Hf+Mwp 



.A 



7 



(A.7) 

All the entries in Hp^ are matrices of infinite rank. The 
single- mode Floquet matrix Hp is on the diagonal and it 
has the famiUar form 



.A 



HM-I/iw hw 

Hl-i] HM + Ifuj 



(A.8) 



7 



In Hp2 (A.7), the fcfiwp-shifted single-mode entries Hp -|- 
Ikhujp are coupled by infinite-rank coupling matrices B^"!, 
defined as 



B[±il 



.A 



Blii^o] 
B[±i^ol 
B[±i'Ol 



(A.9) 



7 



By solving the two-mode Floquet eigenvalue prob- 
lem (A. 5), one obtains the quasienergies and the 



quasienergy states. The energy difference = e+ — e_ 
between two consecutive quasienergies is plotted in 
Fig. 9(a). By applying the periodicity, the single-mode 
quasienergy structure is reconstructed almost everywhere, 
visuaHzed in Fig. 9(b). At the locations where the weak 
probe is in resonance with the single-mode quasienergy 
states (Aq = hwp or Aq = hio — hujp ), a gap, i.e. an 
anti-crossing, opens in between degenerate single-mode 
quasienergy levels, shown in Fig. 9(c). The gap at the 
anti-crossing is the largest wlien it corresponds to a single- 
probe-photon resonance [faint gray lines in Fig. 9(a)]. The 
gaps at the other anti-crossings are opened by increasing 
the probe amplitude, corresponding to the possibility of 
multi-photon probe processes. 

The comparison of the generalized quasienergies [see 
Fig. 9(c)], calculated with Ap/Huj — 0.20 (solid) and 
Ap/Hoj = 0.40 (dashed), gives an example how the probe 
field starts to interplay with tlie single-mode quasienergy 
levels as the probe amplitude ^p increases. By compar- 
ing the two-mode quasienergies calculated with different 
probe ampHtudes, one observes a horizontal shift in the 
location of the anti-crossing, and an enhanced deviation 
from the single-mode quasienergy (dash-dotted). These 
are examples of quantitative deviations from the pertur- 
bative results (13). This kind of a comparison gives a 
qualitative method to study non-perturbatively the higher 
order processes in the probe amplitude Ap. 

The vertical shift of the probe resonances in Fig. 9(c) 
is understood as a Bloch-Siegert ' -type correction due 
to the moderately strong probe field. Moreover, the 
increasing probe amplitude generates effects similar to 
the dynamic (ac) Stark and generalized Blocli-Siegert''' '" 
shifts, but now in terms of the perturbed single-mode 
quasienergy levels. 
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